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Abstract 

In the first part of this article [8], we proved a local version of the circular law up to the finest 
scale N~ 1/2+£ for non-Hermitian random matrices at any point z £ C with \\z\ — 1| > c for any c > 
independent of the size of the matrix. Under the main assumption that the first three moments of the 
matrix elements match those of a standard Gaussian random variable after proper rescaling, we extend 
this result to include the edge case \z\ — 1 = o(l). Without the vanishing third moment assumption, we 
prove that the circular law is valid near the spectral edge \z\ — 1 = o(l) up to scale A r_1 ^ 4+e . 
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1 Introduction 



The circular law in random matrix theory describes the macroscopic limiting spectral measure of normalized 
non-Hermitian matrices with independent entries. Its origin goes beck to the work of Ginibre |16| , who found 
the joint density of the eigenvalues of such Gaussian matrices. More precisely, for an TV x TV matrix with 
independent entries ~~~}^ z ij such that Zij is identically distributed according to the measure /i g = ^e - ' 2 ' dA(z) 
(dA denotes the Lebesgue measure on C), its eigenvalues ixi, . . . , /ijv have a probability density proportional 
to 

n^-^e-^l^l 2 (1.1) 

%<i 

with respect to the Lebesgue measure on C . These random spectral measures define a determinantal point 
process with the explicit kernel (see [16]) 
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with respect to the Lebesgue measure on C. This integrability property allowed Ginibre to derive the circular 
law for the eigenvalues, i.e., j^Pi^ converges to the uniform measure on the unit circle, 



7T 



kl<i 



dA(z). 



(1.3) 



This limiting law also holds for real Gaussian entries |12 , for which a more detailed analysis was performed 
in (7j[l5}[26]. 



For non-Gaussian entries, Girko 17 argued that the macroscopic limiting spectrum is still given by 



( |1.3[ ). His main insight is commonly known as the Hermitization technique, which converts the convergence 
of complex empirical measures into the convergence of logarithmic transforms of a family of Hermitian 
matrices. If we denote the original non-Hermitian matrix by X and the eigenvalues of X by /ij, then for any 
C 2 function F we have the identity 



1 N 1 r 

i vE^) = ^/AF(,)Txlog(^ 

3=1 J 



z*)(X - z)dA(z). 



(1.4) 



Due to the logarithmic singularity at 0, it is clear that the small eigenvalues of the Hermitian matrix 
(X* — z* ) (X — z) play a special role. A key question is to estimate the small eigenvalues of (X* — z* ) (X — z) , 
or in other words, the small singular values of (X — z). This problem was not treated in [17], but the gap was 
remedied in a series of papers. First Bai |3| was able to treat the logarithmic singularity assuming bounded 
density and bounded high moments for the entries of the matrix (see also (4|). Lower bounds on the smallest 



Pan, Zhou (2l] 

Tikhomirov I IS] weakened the moments and smoothness assumptions for the circular law, till 



singular values were given in Rudelson, Vershynin 24[ 25 
and Gotze 



and subsequently Tao, Vu |28 



29 



the optimal L assumption, under which the circular law was proved in 

In the previous article [8], we proved a local version of the circular law, up to the optimal scale 7V~ 1//2+£ , 
in the bulk of the spectrum. More precisely, we considered an N x TV matrix X with independent reaQ 
centered entries with variance TV -1 . Let fij, j £ [[1,TV] denote the eigenvalues of A. To state the local 
circular law, we first define the notion of stochastic domination. 



1 For the sake of notational simplicity we do not consider complex entries in this paper, but the statements and proofs are 
similar. 
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Definition 1.1. Let W be a random variables and ^ be a deterministic parameter. We say that W is stochas- 
tically dominated by if for any a > and D > we have 

\W\ > N a y] ^ N- D 

for sufficiently large N . We denote this stochastic domination property by 

W -< or W=0^(V). 

In this paper, as in (8], we assume that the probability distributions of the matrix elements satisfy the 
following uniform subexponential decay property: 

sup pflVjVJfij-l > A) < ■Or 1 e- x> (1.5) 

for some constant i9 > independent of N. This condition can of course be weakened to an hypothesis 
of boundedness on sufficiently high moments, but the error estimates in the following Theorem would be 
weakened as well. 

Let / : C — > K be a fixed smooth compactly supported function, and f Zo ([i) — N 2a f(N a (fi — z )), where 
zq depends on N and ||zo| — 1| > t for some r > independent of N, and a is a fixed scaling parameter in 
[0, 1/2]. Theorem 2.2 of (8] asserts that the following estimate holds: 



z)J -< N- i+2a . (1.6) 

This implies that the circular law holds after zooming, in the bulk, up to scale N~ 1 / 2+£ . In particular, there 
are neither clusters of eigenvalues nor holes in the spectrum at such scales. 

We aim at understanding the circular law close to the edge of the spectrum, i.e., \zq\ — 1 = o(l). The 
following is our main result. 

Theorem 1.2. Let X be an N x N matrix with independent centered entries of variances 1/N and vanishing 
third moments. Suppose that the distributions of the matrix elements satisfy the subexponential decay property 
(1.5). Let f Zo be defined as previously and D denote the unit disk. Then for any a € [0, 1/2] and any z G C, 
we have 

_ , \ (i7) 

Notice that the main assertion of (1.7 1 is for \zq\ — 1 — o(l) since the other cases were proved in |8 , stated 
m |L6| . 

Without the third moment vanishing assumption, we have the following weaker estimate. This estimate 
does not imply the local law up to the finest scale N~ 1 / 2+e , but it asserts that the circular law near the 
spectral edge holds at least up to scale A r_1 / 4+e . 

Theorem 1.3. Suppose that X is an N x N matrix with independent centered entries, variance l/N, satisfying 
the subexponential decay property (1.5). Let f ZQ be defined as previously, with \zq\ C, and D denote the 
unit disk. Then for any a € [0, 1/4] and any z 6 C, we have 



U-^fM-^ J D f Z0 {z)dk(z)^ -< iV- 1+2Q . 



(1.8) 
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Shortly after the preprint [8] appeared, a version of local circular law was proved by Tao and Vu 30 
under the assumption that the first three moments matching a Gaussian distribution both in the bulk and 
near the edge. Both the assumptions and conclusions of Theorem 20 (30], when restricted to near the edge, 
are thus very similar to Theorem |1.2| On the other hand, in the bulk case the assumption of vanishing third 
moments were not needed in [8] . Our proof in this paper follows the approach of the companion article (8] , 
except that for small eigenvalues we will use the Green function comparison theorem |13| . On the other 
hand, the method in |30| relies on Jensen's formula for determinants. 

A main tool in the proof of the local circular law in [8] was a detailed analysis of the self-consistent 
equations of the Green functions 

G ij (w) = [(X*-z*)(X-z)-w]r j 1 . 

We were able to control Gij(E + irj) for the energy parameter E in any compact set and for sufficiently 
small r\ so as to use the formula ( |1.4[ ) for functions F at scales A r_1 / 2+e . We proved that, in particular, the 
Stieltjes transform m = iV _1 Ga converges to m c (w, z), a fixed point of the self-consistent equation 

l + wm(l + m) 2 

g( m ) = r— la — i (1-9) 

\z\ 2 - 1 

when ||z| — 1| ^ r > 0. In this region of z, the fixed point equation is stable. However, for \z\ — 1 = o(l), 



(1.9) becomes unstable when \w\ is small. In this paper, we will show that m is still close to the fix point 
of this equation under the additional condition that the third moments of the matrix entries vanish. With 
this condition, we can compare m of our model with the corresponding quantity for the Gaussian ensemble. 
In order to carry out this comparison, we prove a local circular law for the Ginibre ensemble near the edge 



\z\ — 1 = o(l) in Section 4. We now outline the proof of the main result, Theorem 1.2 Step 1: we use the 
Hermitization technique to convert the Theorem 1 1 ,2| into the problem on the distribution of the eigenvalues 
Aj's of (X — z)*(X — z). Step 2: for the eigenvalues of order 1, i.e. A, ^ e for some e > 0, we will control 
^log(Ai) via an analysis of the self-consistent equations of the Green functions 

G ij {w) = [{X*-z*){X-z)-w]T l . 

Step 3: for the eigenvalues of order o(l), we will show that ^log(Ai) of our model has the same asymptotic 
distribution as the Ginibre ensemble via the Green function comparison method first used in |13| . Step 4: 
combining the previous steps, we show a local circular law for the Ginibre ensemble. 



2 Proof of Theorem [L2] and Theorem 11.31 

2.1 Hermitization. In the following, we will use the notation 

Y Z =X - zl 

where / is the identity operator. Let Xj(z) be the j-th eigenvalue (in the increasing ordering) of Y*Y Z . We 
will generally omit the z— dependence in these notations. Thanks to the Hermitization technique of Girko 



17 , the first step in proving the local circular law is to understand the local statistics of eigenvalues of Y*Y Z . 



More precisely, (1.4 1 yields 

N-^fM = ^ K N ~ 1+2a J (A/)(£)5^1ogAj(.z)cIA(£), (2.1) 
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where z — zq + N~ a £. Roughly speaking, for any fixed e > 0, the sum iV -1 J2\ > £ \°g^i( z ) wm approach 
J x>£ (\ogx)p c (x)dx, where p c is the limiting spectral density of Y*Y Z . To estimate the error term in this 
convergence, we need to know some properties of p c (Subsection 2.2) and, most importantly, the strong 
local Green function estimates of Subsection 2.3 The \ ^ e part will be studied with the Green function 
comparison method in Section 3. 

2.2 Properties of m c and p c . Define the Green function of Y*Y Z and its trace by 

1 1^1 

G(w) := G(w, z) = (Y* Y z — w) , m(w) :— m(w, z) = — Tr G(w, z) = — } , w = E + irj. 

W TV ^ — ' Xi(z) — w 

j= i ^ > 

We will also need the following version of the Green function later on: 

G(w) :=Q(w,z) = {Y z Y;-w)-\ 

As we will see, for some regions of (w, z). with high probability m(w, z) converges to m c (w, z) pointwise, as 
N — >• oo where m c (w, z) is the unique solution of 

to" 1 = -w(l + m c ) + |z| 2 (l + m c )" 1 (2.2) 

with positive imaginary part (see Section 3 in |18| for the existence and uniqueness of such a solution) . The 
limit m c (w, z) is the Stieltjes transform of a density p c (x, z) and we have 

m c ( w ,z) = [ ^^dx (2.3) 



whenever r\ > 0. The function p c (x,z) is the limiting eigenvalue density of the matrix Y*Y Z (cf. Lemmas 
4.2 and 4.3 in (3j). This measure is compactly supported and satisfies the following bounds. Let 

A±:=A±W:= 8^TlT' a: =v/l+W- (2-4) 

Note that A_ has the same sign as \z\ — 1. It is well-known that the density p c (x,z) can be obtained from 
its Stieltjes transform m c (x + vq, z) via 

Pc(x,z) = -Im lim m c (x + it], z) = -l x e[max{o,A_},A + l Im lim m c (x + iry, z). 

TT j;-S-0+ 7T IJ-S.0+ 

The propositions and lemmas in this subsection summarize the properties of p c and m c that we will need in 
this paper. They contain slight extensions (to the case \z\ = 1 + o(l)) of the analogous statements 18]. We 
omit the proof, strictly similar to the calculation mentioned in j8j. 

In the following, we use the notation A ~ B when cB ^ A ^ c -1 -B, where c > is independent of N. 



5 




Proposition 2.1. The limiting 
density p c is compactly sup- 
ported and the following prop- 
erties hold. 

(i) The support of p c {x,z) is 
[max{Q,A_},A+]. 

(ii) As x — > A + from below, 
the behavior of p c {x,z) 
is given by p c (x,z) ~ 

(Hi) For any fixed e > 0, if 
max{0,A_} + e ^ x < 
A+ — e, then p c (x, z) ~ 1. 

(iv) Near max{0, A_}, if z is 
allowed to change with 
N such that \z\ = 1 + 

En (with en — o(l)J, Figure 1: The limiting eigenvalues density p c (x, z) for z = 1.5, 1.3, 1.2, 1, 0.4 
i/ien A_ = o(l) and i/ie 
behavior of p c (x,z) de- 
pends on En- We will not need this detailed property and will only need the following upper bound: 
for any 5 > there is C > such that for any \z\ S^ 1 , x G [0, 5], we have p c (x, z) ^ C/y/x. 



From this proposition, we obtain the following bound regarding the Stieltjes transform of p c (see (2.3 1). 
Proposition 2.2. Uniformly in z and w in any compact set, the Stieltjes transform of p c is bounded by 

m c {w,z) = 0(M~ 1/2 ). (2-5) 



Moreover, for any fixed £ > 0, uniformly in £ < \w\ < £ , \z\ < £ 1 , we have m c 



1. 



We now collect some technical properties of m c used in this paper. Define k := k(w, z) as the distance 
from E to {A + , A_}: 

K = vam{\E-X-\,\E-X + \}. (2.6) 

For \z\ ^ 1, we have A_ < 0, so in this case we define k := \E — A+|. 

The following two lemmas, concerning the case c ^ \w\ c _1 , are analogous to Lemma 4.1-4.3 of [Si, 
and the proofs are essentially the same. Notice that the properties of p c {x) used in 18] when c ^ \w\ ^ c _1 
can be summarized as follows: (1) p c {x) ~ \/A+ — % if < c < x < A+. (2) J° (t] p c (x) = o(l). From 
Proposition 2.1 these two estimates hold uniformly for all \z\ ^ C, including \z\ — 1 = o(l). We therefore 
can use the proofs in (8] to obtain the following two Lemmas. 

Lemma 2.3. There exists tq > such that for any t SC to, if \\z\ — 1| ^ r and r SC |ui| t _1 £/ien i/ie 
following properties concerning m c hold (all constants in the following estimates depending on t). Notice 
that by fl^ , we have in this case A_ = fej^ = 0(||z| - 1| 3 ) = 0(t 3 ). 



G 



Case 1: E ^ A + and \w — A+| ^ r. We have 

| Re m c | <- 

Case 2: \w — A+| ^ r. We have 

9 

m c (w,z) 

and 



1, ^ Rem c ^ , Imm c ~ n. 



3 + a 



/8(l + a) 3 
a(3 + a) 5 



- A+) 1 / 2 + 0(A + - w), 



(2.7) 



(2- 



Inn 



^= if k ^ rj and E ^ A+, 
^ if k ^ rj or E ^ A+. 



Case 5: |u; — A+| ^ r and i? A + (notice that E is allowed to be smaller than \_). We have 

|m c | ~ 1, Imm c ~ 1. 



(2.9) 



(2.10) 



Lemma 2.4. There exists tq > smc/i that for any t ^ To if the conditions \\z\ — 1| ^ r and t ^ \w\ ^ r _1 
ZioZd, f/ien we have the following three bounds concerning m c (all constants in the following estimates depend 
on t ): 

|m c + l| ~ \m c \ ~ M~ 1/2 , (2.11) 



Im 



(-i + |z 2 |) k- 



i 



w(l + m c ) 
-2 
3 + a 



si CImm c , 

-2 

i c 

3 — a 



w 



(2.12) 
(2.13) 



2.3 The Strong Green function estimates. We now state precisely the estimate regarding the convergence 
of m to m c . Since the matrix Y*Y Z is symmetric, we will follow the approach of |14| . We will use extensively 
the following definition of high probability events. 



Definition 2.5 (High probability events). Define 

ip := (logA0 loglogAr . 



(2.14) 



Let £ > 0. We say that an N -dependent event Q holds with ^-high probability if there is some constant C 
such that 

P(fi c ) < 7V c exp(-/) 

for large enough N . 

For a ^ 0, define the z-dependent set 

S(a) := {weC : max(A_/5,0) <^<5A + , ^iV-Vcr 1 «S»7< 10}, 
where ip is defined in ( 2.14[ ). Here we have suppressed the explicit ^-dependence. 
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Theorem 2.6 (Strong local Green function estimates). Let e > be given. Suppose that z is bounded in C, 
uniformly in N. Then for any £ > 0, there exists > such that the following event holds with (-high 
probability: 

p| h m (w)-m c (w)\^^ c <^\. (2.15) 

weS(C ( ),\w\>e 

Moreover, the individual matrix elements of the Green function satisfy, with (-high probability, 



H { V \G v -m oSv \< V <><U^ + jL)\. (2.101 



'j 



weS(C ( ),\w\>e 

Proof. The proof mimics the one of Theorem 3.4 in |8], in the case of general z but restricted to \w\ > e. 
The proof is the same as under the restricted hypotheses \ \z\ — 1| > r, |iy| ^ e, a special case satisfying the 
hypotheses of Theorem 3.4 in (8]: in our proof for this special case, we used the fact that r is a positive 
constant very indirectly. We only required some properties of m c in this case, which are Lemmas 4.1 and 
4.2 in |8J, and have as analogues Lemmas 2.3 and 2.4 from the previous subsection. □ 

2.4 Conclusion. Following the Hermitization explained in the previous section, in order to prove Theorem 



1.2 we need to properly bound (2.1). Now logAj(z) needs to be evaluated even for z close to the edge, and 
for this we first define a proper approximation of the logarithm, with a cutoff at scale N~ 2+e . 

Definition 2.7. Let h(x) be a smooth increasing function supported on [1, +oo] with h(x) ~ 1 for x ^ 2 and 

h(x) = for x ^ 1. For any e > 0, define <f> on R + by 

<P(x) = <j> £ {x) = h(N 2 - 2e x) (logx) U-h (^)) • (2-17) 

We now prove that is a proper approximation for log in the sense that 

Tr (log(F*y) - <p(Y*Y)) -< N Cs (2.18) 

uniformly in z in any compact set, where C is an absolute constant independent of N. 

For this purpose, we first bound the number of small eigenvalues of Y*Y. This is the aim of the following 
lemma, proved in Subsection |3.2| 



Lemma 2.8. Suppose that \w\ + \z\ ^ M for some constant M independent of N . Under the same assumptions 
as Theorem 1.3 for any ( > 0, there exists Cq such that if 

r] ^ ip c <N- l \w\ 1/2 (2.19) 

then we have 

\m(w,z)\ sC (logA0|wr 1/2 (2-20) 
with (-high probability (notice that if we take rj ~ \w\ then the restriction is \w\ 
Using 

\{j : Xj G [E-r),E + rj\}\ < Cr]lm.m(E + ir]) (2.21) 
and Lemma 2.8 with E = and 77 = N~ 2+2s , we obtain 

1|{j:|A,K7V- 2 + 2£ }|^. (2 22) 

Moreover, we have the following lower bound on the smallest eigenvalue. 



<S 



Lemma 2.9. Under the same assumptions as Theorem \1.3\ 

| log Ax (0)| -< 1 

holds uniformly for z in any fixed compact set. 



Proof. This lemma follows from 25 or Theorem 2.1 of |28], which gives the required estimate uniformly 
in z. Note that the typical size of Ai is N~ 2 25 , and we need a much weaker bound of type P(Ai(z) ^ 



') ^ N c for any e, C > 0. This estimate is very simple to prove if, for example, the entries of X have 



a density bounded by N c , which was discussed in Lemma 5.2 of 



The preceding lemma together with (2.221 implies 

Tr 



(log(y*Y) - 4>(Y*Y)\ -4 N Ce , 



□ 



(2.23) 



where <fi(x) = h(N 2 2e x)(logx). The proof of (2.181 is then complete thanks to the following easy lemma 
bounding the contribution of exceptionally large eigenvalues. 



Lemma 2.10. Under the same assumptions as Theorem \1.3\ for any e > 0, 

Tr U(Y*Y) - <j>(Y*Y)\ -< 1 
holds uniformly for z in any fixed compact set. 



(2.24) 



Proof. Notice first that, for any cq > 0, there is a constant C > such that, uniformly in \z\ < cq. Xn ~< N c 
(this is for example an elementary consequence of the inequality An ^ Tr(Y*Y)). The proof of the lemma 
will therefore be complete if 

\{i : Xi > 2A+}| -< 1. (2.25) 

Let {X^)keli,N} De the increasing eig envalues of y(*)*yW, where Y^ = X - (1 - S)z. By the Weyl 
inequality for the singular values, for any 5 > 



Af^Ai 5)1/2 + ( 5. 



(2.26) 



Moreover, as |(1 — S)z\ < 1, we are in the context of Lemma 5.1 in |8j, which yields, for any £ > 0, the 
existence of > such that 



A 



("5) 

'N-(p^< 



^ A 



(6) 



(2.27) 



with C-high probability. Equations (2.261 and (2.271 give 



X N-<p°i ^ A 



(«5) 



26\ { l )1/2 + 5 2 . 



Together with A^? — >• A + as 5 — >• (uniformly in \z\ < cq), this concludes the proof of (2.251 by choosing 5 



small enough. 



□ 



Remark 2.11. With some more technical efforts, the techniques used in this paper allow to prove the stronger 
result Ajv — A + -< N~ 2 ^ 3 . This rigidity on the edge could be proved similarly to the rigidity of the largest 
eigenvalues of X* X , for non-square matrices |22|. 
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It is known, by Lemma 4.4 of j3], that 



(logx)A z p c (a;, z)dx = 4xd(z) 



By definition of 4> an d p c 



x)p c (x)dx - / \og(x)p c (x)dx 



< N~ 1+Cs 



Then from equations (2.1 1 and (2.18), Theorem 1.2 follows if we can prove estimate 



A/(£) ( Tr <b(Y*Y)-N / 4>{x) p c {x)dx ) d£d£ 



(2.28) 



(2.29) 



(2.30) 



We first evaluate Tr <fi(Y*Y) in the following lemma. 



Lemma 2.12. Let x be a smooth cutoff function supported in [—1,1] with bounded derivatives and x{v) = 1 
for \y\ 1/2. For the same e presented in the definition of <j) e , we define the domain in C 



I = iN- l+e \/E < n, \w\ ^ s\ 



There is some C > such that for any e > 0, (recall: (j) := <j) e from Def. 2.7) 

Ti(f)(Y*Y)- N j <j)(x)p c (x)dx - — ( x{v)<f>'(E)Re(m{w) -m c {w))dEdT] 



-< N' 



Cc 



(2.31) 



(2.32) 



uniformly in z in any compact set. 



Proof. Suppose that (2.321 holds with the integration range / replaced by / where 

1 = {N- 1+S VE sC r)\ . 

Define the notation Am = m — m c . One can easily check that if w € / and (f>{E) ^ 0, then n ^ N~ 1+e \w\ 1 ' 2 . 
Using Proposition |2.2[ we have 



i/i 



x(ri)4>'(E) Re (Am(iu)) dE<drj < 



e<|w|<c,w- 1 + c Vb<77<c NEr/ 



Hence (2.321 holds with integration range given by /. 



We now prove (2.321 with / replaced by /. The Helffer-Sjostrand functional calculus (see e.g. 

iyq"(x)x{y) + Kq( x ) + W \ x ))x' \y) 



for any smooth function q : 



11 



gives 



3(A) = 



2tt 



A 



x-iy 



-dxdy. 



As a consequence, 



Tr cj>(Y*Y) = — / (ir](j)"(E)x(v) + i</>{E) X '{v) - #'(£)x'(l))m(£ + i^dEdr), (2. 

* Jr,>Q V > 



33) 



N / <j){x) p c {x)dx 



N 



vn<t>"{E) X {ri) + ^(E)x'(v) ~ V<l>'(E)x'(ri))m c (E + i^dEdr,. (2.34) 
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By definition, x'iv) 7^ on ly if \v\ € [1/2, 1]. Using (2.151, we have 
N 



4(E) X '(rj) - rj^'(E)x'(v) Am(£ + ir])dEdr] -< 1. 



hence the difference between (2.331 and (2.34) is 

f N 

Tr 4>{Y*Y) - TV / (f){x)p c (x)dx = — 



ir)<p"(E)x(r))Am(E + ir/)dE dr? + (1). 



r;>0 



Tr</>(F*Y) -iV / 0(a;)p c (a;)da: 



r](j)" (E)x{v) lmAm(E + \r])dEdr] + 0^(1) 



Since the left side of the equation, x an d <fi are a U real, this can be rewritten as 

N 

Furthermore, for r\ ^5 N~ 1+£ \/E and E ^ 7V~ 2+2e , we have i]Imm(E + irj) ~< TV -1 , because 77 1— > 
rjlmm(E + irf) is increasing, and at 770 = N~ 1+6 y/E the result its true thanks to ( |2.20[ ) (77 ^ N~ 1+£ \fE 



and E ^ N~ 2+2e imply that rj > 



w 



1/2 



b"(E)\ < C(l + I log(S)|)£:- 2 , the above estimate yields 
iV / 7 ? 0"(£;)x(r7)Imm(w)d J Bd77 ^ 



as we saw previously). Together with the easy bound 



CN -l+Ce E -3/2 dE ^ N Ce^ 



/AT- 2 ^_E^O(l) 

The same inequality holds when replacing m with m ci hence 



(2.35) 



N rjcj)" (E)x{r{)lm{Am{w))dEdri < N 



Ce 



To estimate N fr ?]<f>" {E)x(v) hn{Am(w))dEdr], we integrate this term by parts first in E. then in 77 (and 
use the Cauchy-Riemann equation -J^ Im(Am) = — Re(Am)) so that 



N I jqcj)" (E)x{ri)lm{Am{w))dEdri =N J r/x(r])(j)' [N 2 ~ 2e r] 2 )l-ai{Am(N 2 - 2£ r] 2 + ir/))dT] 

- N J x{v)<P'i E ) Re(Am(7i;))d£;d7; 

— N [ T]x'(vW(E)'Rie(Am(w))dEdr]. (2.36) 



The last term is nonzero only if |?y| € [1/2, 1]. By (2.151, this term is of order 0^(1). 

Concerning the first two terms, they can be bounded in the follow ing way: as we already saw, if w € I 
and <fi(E) ^ 0, then 7/ ^ N- 1+E \w\ 1/2 ; one can therefore use \2.20\ : \m\ + \m c \ < (log TV) |w | _ 1/2 . This 
together with \<fi'(E)\ ^ C(l + | \og(E)\)E^ 1 proves that the first two terms are 0^(N Ce ), completing the 
proof: only the third term remains. □ 



To complete the proof of Theorem 1.2 we estimate the third term in (2.36 1 under the additional assump- 
tion that the third moments of the matrix entries vanish. The following lemma provides such an estimate 
whose proof will be postponed to Section [373] To state this lemma, we introduce the following notation, 

.(9) _ 



= N / Ag(0 / X (v)^(E)Re(m 1 (w) - m 2 (w))dEd V dA^), 
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(1)* (1) 

for any given smooth function g and initial random matrices X\ , X2 , with associated matrices Y£ Y z , 

(2)* (2) 

Yz Y z having respective Stieltjes transforms mi and rri2- Moreover, we will write 



7 (a) 

J X,c 



N I Ag(£) J x{vW(E)Bje(m(w) - m c {w))dEdr)dA(£). 



Notice that x ^ and z| c depend on e through the definition of 4> in Def. 



2.7 



Lemma 2.13. Under the assumptions of Theorem \l.S\ there exists a constant C > such that for any small 
enough e > 0, we have 



7 U) 
J X,c 



where Cf is a constant depending only on the function f . 



Combining this lemma with (2.321, we obtain (2.301 and complete the proof of Theorem 1.2 



To prove Theorem |1.3| instead of Lemma |2.13| we only need to prove the following lemma which does 
not assume the vanishing third moment condition. This lemma will be proved at the end of Section 3. This 



concludes the proof of Theorem 1.3 



Lemma 2.14. Under the assumptions of Theorem \1.3\ for some fixed C > 0, for any small enough e > 0, we 
have 

Z [ H ■< N^c f , 
where Cf is a constant depending only on the function f. 



3 Proof of Lemmas 12.81 and 12.131 



3.1 Preliminary lemmas. This subsection summarizes some elementary results from [8], based on large 
deviation estimates. Note that all the inequalities in this subsection hold uniformly for bounded z, no 
matter its distance to the unit circle. We first introduce some notations. 

Definition 3.1. Let T,U C [1,AT]. Then we define Y^^ as the (N - |U|) X (N - |T|) matrix obtained by 
removing all columns of Y indexed by i £ T and all rows of Y indexed by i € U. Notice that we keep the 
labels of indices ofY when defining y( T < u ). 

Let yi be the i-th column ofY and y t be the vector obtained by removing yi(j) for all j £ S. Similarly 
we define y^ be the i-th row ofY. Define 



£(T,U) 



(y(T,Uh* F (T,U) _ w 
y(T,UW(T,lim* _ w 



(T,U) 



(T,U) 



TrG (T,U) 



1 

~ N 

= l T r£( T < u ). 
N 



By definition, m@'^ = m. Since the eigenvalues ofY*Y and YY* are the same except the zero eigenvalue, 
it is easy to check that 



(T,U) 



(w) — mg + 



(T ,u) , |U|-|T| 



Nw 



For |U| = |T|, we define 



m (T,U) . = m (T,U) =m (T,U) 



(3.1) 
(3.2) 
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Lemma 3.2 (Relation between G, G^'^ and G^' T ^). For i,j^k ( i = j is allowed) we h 



ave 



GikGkj M$,k) n _ QikQkj 



and 



r m r , (Gy*)(y l G) (0>i) _ (G^y*) (y»G^>) 

+ 1-yiGyj ' G ~ G " i + y,G(M y: ' 

0(i ,B) _ c , (gyQ (y?g) c _ e( i,0) _ (g (i - 0) yi)(yi*g (a) ) 



(3.3) 
(3.4) 



Furthermore, the following crude bound on the difference between m and m G ' ' holds: for U,Tc [1, NJ 
we have 

\m - mg™| + \m- T) | < |U| + |T| . (3.5) 
y Nn 

Definition 3.3. In the following, Ex means the integration with respect to the random variable X. For any 
T C [1, NJ, we introduce the notations 



and 



Z™ := (1-EyJyf^yp 
4 T) Ml-EyJyP^V^. 



Recall by our convention that y^ is a N x 1 column vector and yj is a 1 x TV row vector. For simplicity we 
will write 

z i = zf\ z i = zf\ 

Lemma 3.4 (Identities for G, Q, Z and Z). For any Tc [1, AT], we have 



1 + ™^ + !^+^ 



G 



(0,T) 
ij 



where, by definition, Q^'"^ = if i € T. Similar results hold for Q: 

-l 



(T,0) 



C (T,0) 



= —w 



(3.6) 
(3.7) 

(3.8) 
(3.9) 



Lemma 3.5 (Large deviation estimate). For any ( > 0, there exists > such that for T C [[1,-/V], 
|T| ^ N/2 the following estimates hold with (-high probability: 



,(T), 



\z 



(T), 



(l-!E yi )(yf)G^)yp) 



Immg' 1 ' + \z\ 


»ImG<™ 


Nn 


jlm.mg' 1 " > + \z 





(3.10) 



Nr) 
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Furthermore, for i ^ j , we have 

(l-E yiy j(yf>GW%f>*) 



(l-E yiy J yf^« 



,(T)* n (y,T)„(T) 



sC<^2</ 2 1 









N77 


Immg J ' T ' + 


NI 2 Im^' T) 4 


M 2 Imejf T) 



Nr] 



(3.11) 
(3.12) 



where 



E 



(yf^Vf)*) = 1,^), Eyiyj (yf)^W.t) y ro) = 1,1^. (3.13) 



3.2 Proof of Lemma 2.8 By (|3.10|) and (|3.5|, with (J-high probability, we have for some large constant C 

1 



Nr] r Nr) Y 



j Imrn| 
Nr] 



1 



Imml 



Nri (log AO c ' 



(3.14) 



where we have used the Cauchy-Schwarz inequality in the last step. Moreover, choosing T = {i} in (3.8), we 
have 

-l 



,(i,0) 



If 



3C 



Afy (log7V) c 



and Iml ^ 3, 



then by (3.151 and (3.51 we have 



,(i,0) 



> hi 



1 i CM) t 

1 + TUq + Z\ 



2 , 



(log TV) c 



(3.15) 
(3.16) 

(3.17) 



Similarly, with £-high probability, we can bound Zi by 



Imml' 1 '' + |z| 2 Im (/,.,*' ^ 



N77 



< tp c < — + 



1 Im to + I z 1 2 Im 0, 



iV?y (log AO c 



(3.18) 



If ( |2.20| is violated, i.e., |m| ^ (log N)\w\~ 1/2 , then ( |3.16| is implied by the assumption fl2.19| . Hence ( |3.17 l 
holds. By ( pTej l, 

I^G.r 1 = |1 + m^ 0) + |z| 2 ^ 0) + Z,| ^ |m|/2 — C— \w\~^ 

and we conclude that \Gu\ ^ | -uj | 1 / 2 with £-high probability, for any 1 sC i ^ AT. This contradicts the 
assumption \m\ ^ (log N^w^ 1 / 2 and we have thus proved (2.201. 

3.3 Proof of Lemma |2.13| We begin with the following estimates on the elements of the Green function, 
weaker than those in Theorem |2.6| but valid for any w and z. 
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Lemma 3.6. Suppose \w\ + \z\ $J C. For any C > 0, there exists Cq such that if the assumption (2.19) holds 
then the following estimates hold: 

max|G 4l | ^(logAOM- 1 / 2 , (3-19) 



maxHIGwH^l^aogJV) 



max | G i:j | C (log Ny\w\ 



2u„|-i/2 



(3.20) 
(3.21) 



with Q-high probability. Furthermore, for eigenvectors u a ofY*Y are delocalized: with (,-high probability 



max|K||^ ^ tpCtN- 1 . 

a 

Proof. From equations (6.37) and (6.38) of [8], the following estimate holdfj^j 



>(i,0) 



w(l + m) 



where 



£ i = w 



* = 



(1 + m) 2 



(M) i ry(i) 

m G m + &i 







\z 



(0|2 



|to||l + m| 3 



= 0(</^ /2 *) 



(3.22) 



(3.23) 



(3.24) 



' Im m c + I m — m c 



Nrj 



Nn 



By using (2.15| > -||2.16| | , fl3.23| ), ( |3.24[ ) and ( |2.1l| , for |w | > g one ca n be easily prove: |Gij| + |S^' 0) | sC 
which implies ( |3.19[ )-(|3.21^ Moreov er, if \ z\ ^ 1 /2, ( |2.15[ )-( |2.16| still hold without the restriction \w\ ^ e 
(see Theorem 3.4 in |8j). It implies ( 3.19 )-( 3.21 1 in the case \z\ ^ 1/2. From now on, we therefore assume 
that for some small enough e > 0, 

M<£, |*| > 1/2. 



-1/2 



From (3.6 1, (3.5) and (3.8), we have 



Nn 



+ W ij +w\Z i \+w\Z j \, 



where 



Wi.. 



\z\ 2 ((1 + m£ l) + Z®)- 1 - (1 + mg j) + Z^)- 1 ) 



(3.25) 



If e 1 / 4 ^!" 1 / 2 < |m| < (logA0M~ 1/2 then by ( |3.14| we have \Z^\ < (logiV)" c |TO| due to the restriction of 
77 in (2.191. Thus we can estimate Wij by 

\W id \ < |z| 2 [j^>| + + (ATt?)" 1 ] (l + mM+Z®)-\l + mU'f> + zV ) )- 1 ( 3 - 2 6) 
^C|*| 2 |l + m|- 2 [izfl + l^l + ^r?)- 1 ! < (logAT)- ^! 1 / 2 (3.27) 



2 Although these bounds in [S] are stated under the assumption \\z\ — 1| ^ r, the same argument holds also for 2 close to the 
unit circle, under the extra assumption \w\ > e. 
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Under the assumption is small enough and e 1 /^]- 1 / 2 < |m| (log N^w]- 1 / 2 , the condition (|3.16| 



is satisfied and, from (3.181, we have 

\Zi\ < 

Thus we have in this case the estimate 



Ntj (logN) c 



-1/2 



\G U - m\ < iV- 1 l G « - %l < l G « G «l l G ^ - G 7/l < l G " G «l 



Nrj (log TV) c 



(3.28) 



|w| 1/2 (3.29) 



Define the parameter 



7(2, u>) = max |G^||w;| 1 ' /2 

i 

Then ( |3.29| ) and the assumption e 1/4 \w\~ 1/2 < |m| < (log A0M~ 1/2 imply that 

j(z,w) 2 



7(2, w) < Glog7V + 



(logiV) c, 



By continuity (in 77) method, we have proved 7(2:, w) ^ Clog TV (i.e., (3.19)) assuming that e 1 / 4 ^! 1 / 2 ^ 
| m| < (log A0M -1/2 holds. Since |m| (logA0M~ 1/2 nolds with C-high probability by ( |2.20[ ), to prove 
(3.191, we only have to consider the last case \m\ ^ e 1 / 4 ^! -1 / 2 . 
By ( |3.15[ ) and ( |3.14| ), we have 

|^' 0) | ^ l^l^ll + m^+Z^r 1 ^ e-^l^wV 1 ! 2 12 

Thus 

i4 l,0) + m 2 ^ ,0) i > e- l/4 kr l/2 /2 - e i/4 kr i/2 - i > e-^lwr^/A 



By ( |3.18[ ) and the notation m^' 0) + |z| 2 ^' 0) = A, + \B t , A,, B, e E, we have, for ?y satisfies ( |2.19| ), that 

(logiV) c 



\Zi\ ^ (logJV)- 1 ^^ 



Nrj 



From (3.6), we have 
1 



\G, 



^ \w\ [\1 + Ai + iBA - \Zi\] > w 



\l + Ai + iBt\ 



(logiV)" 



(3.30) 



^e^'^w^/B (3.31) 



We have thus proved (3.191. 



We now prove (3.20). From (3.6 1, we have 

( W G l4 r 1 + i+-f ) = -[N 2 ^;' 0) +^] 

By the triangle inequality, we have 

K 2 G^ 0) | < | W G 4I [|z| 2 <^ 0) +Zi]\ + \wGnZi 



(3.32) 



(3.33) 
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Together with (3.321 and using (3.181 to bound Zi, with (2.191 we have 



2 <W 0) I < \wG u [1 + m^ } }\ + 1 + <p c < \wG u 1 1 



< \wG H [l + m^]\ + 1 + + <P-^\ W Z*G U G^\ 



' Imm 



l*l 2 Im<# fl) 



(3.34) 



Then with ( |2.20| and (|3_19| on m and Gu, we obtain ( |3.20 l. 
For any eigenfunction with eigenvalues |A Q — E\ ^ r/, we have 

il\u a (i)\ 2 



C llmGiil < \G r 



and thus 



|u a (i)| 2 <27 ? |G«| < ^cjv- 1 , 



where we used (3.191 for r\ = (p° c N~ 1 \w\ 1 ' 2 . Since this bound holds for all energy E, we have proved (3.221 
(notice that the proof of (3.221 for w ^ e follows by the same argument). 
We now prove (3.21). By the Cauchy-Schwarz inequality, we have 



v |A a -i?lK(j)| 2 \\ a -E\\u a (j)\ 2 



(3.35) 



The imaginary parts, ImG^ + ImGjj, can be bounded by ( |3.19 1. The first term on the right hand side is 
bounded by 



\- \\ a ~ E\\u a (i)\ 2 ^ x ^ (\ a -E)\u a (i)\ 2 



(3.36) 



The real part, ReGu is bounded again by (3.191. Let E k = E + (2 fe - 1)77 and rj k = 2 k 1 rj. From (3.191, we 
have 

CT] k \u a (i)\ 2 
{\ a -E k y + ril 



v {\ a -E)\u a (i)\* C ^ N v 
^ (\ a -E) 2 +ri 2 " ^ ^ 



, 2 2 ^ciogiv^r 1 / 2 



(3.37) 



This concludes the proof of (3.21 ) and Lemma 3.6 



□ 



Lemma 3.7. Suppose that \w\ + \z\ = 0(1). For any ( > 0, there exists Cq such that if (2.191 holds then we 
have 

(3.38) 



and 



and 



and 



\[YG(z)} kl \ + \[GY*(z)} kl \^C(\ogN) 2 
\[YGY*{z)] kl -5 kl \^C{\ogN?\w\V* 



max 

ij 



[G 2 ]ij\ < max 1 1 G\% I < 



C(logAT) 2 |tti _ 1/2 



max j [FG 2 ] y | < 



C(logAT)- 



(3.39) 
(3.40) 

(3.41) 
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and 



with Q-high probability. 



max [FGTjij < 



C{logNf Hl/2 



(3.42) 



Proof. We first prove (3.38) and (3.391. From the general identity 



A(A*A - w^A* = 1 + w{AA* - w)' 1 , 



we have 



(YGY*)ij = 5ij + wQij 
This proves (3.391 with (3.191. To prove ( |3.38 I, we notice first the identity 

YG = QY, GY* = Y*g. 

Recall the rank one perturbation formula 

_ x (A^v^ivA- 1 ) 



and Lemma [3J2J We have 



Together with the equation 



(A + vV)- 1 = A 



Gyi 



1 + vA-'v* 



(see Equation (6.18) in fSl for a derivation) and the definition of y; as the i-th column of Y, we have 



{Qyi)i = -wGuig^y^ = -wzGnQ 3i 
From the large deviation lemma, we can bound 



{%,$) 



(3.43) 



N-q n Nrj 



with high probability. Hence 



\{gY) 3 i\ = \{Gy t )j\ < l^lwzGug^l + CilogN) 4 



(3.44) 



(3.45) 



where we have used (3.191, (3.201 and (2.191. Now we estimate Qjl , with (3.6l, (3.11|, (3.13), we have 



wG u g 



w 2 G lt g^g^ 



»M "33 



! Imm% ij) + \z\ 2 ImGff j) 



(3.46) 



Then using wG»g t ( -' 0) < (log N) 4 and wG)*'*' Gff J > (logiV) 4 (see p^Ofr ), we obtain it is less than (logiV) 2 
This proves KM. 



n (i,ij) 
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By the Cauchy-Schwarz inequality, we have 



|Gl!<./£l G d 2 wEl G ^l 2 < 



' Im Gjj / Im G 



J 3 



V 



V 



which implies (3.401. Using \w\ ^ Imw — rj, (3.211 and (3.40), we have 

\{yg 2 y) 13 \ = \g ij + w gU^c\w\- 1 ' 2 + c\w\- 1 ' 2 r 1 - 1 ^cm^'V 1 



This proves (3.421. Notice that we claimed that all estimates proved for G are valid for Q as well. This is 
because that we can write YY* in Q as A* A where A — Y* . Hence all estimates hold for G will hold for Q. 
By the the Cauchy-Schwarz inequality, we have 



\{YG%A < \Yj?G)ikG ki \ < \(YGG*Y*) ii \\(G*G) i 



Together (3.40) and (3.421, we have proved (3.411. 



□ 



Remark 3.8. The previous lemmas 3.6 and 3.1 also hold if one entry of X is supposed to vanish (this will be 
useful for us in the following of this subsection). Indeed, we just used the following facts: the independence of 
the entries, the subexpotential decay, EAT a 6 = 0, and E|AT a h| 2 = 1/N. We notice that only this last condition 
is changed, and it was are only (3.6 I and ( |3.8| . Furthermore, for ( |3.6[ ), this affects the case i = a as 

-l 



G (0,T) 



1 + m 



(a,T) 



1 



N @bb 



(o,T) 



The difference is therefore just j^Q^' T ^ which is of order 1/N of rn~g' T \ so the proof holds for such matrices 
with a vanishing entry. 

The following lemma is not useful for the proof of Theorem |2.13| but it helps to understand how the 
Green's function comparison method works at here. Furthermore, some by-products are very useful for the 
whole proof of Theorem |2.13| 



Lemma 3.9. Suppose that \w\ + \z\ < c for some fixed c > 0, and that r/ satisfies (2.19). Assume that we 
have two ensembles X\,Xi, both of them satisfying (1.5), and with matrix elements moments matching up 
to order 3. Then we have 



|E X (i) (m(w, z)) - E x( 2)(m(w, z))| 



Nrj' 



(3.47) 



Proof. For k e JO, A^ 2 ], define the following matrix X k interpolating between X*- 1 ' and X^: 

X k (i,j)-- 

Note that X^ — Xq and X^ = Xj^2. A sufficient condition for (3.471 is that for any k ^ 1 



X^(i,j) if k<N(i-l)+j 
X<V(i,j) if k^N(i-l)+j 



\Ex k m(w, z) — ¥,x k _ 1 m(w, z) I ^ 



<P 



N 3 rj 



(3.48) 
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We are going to compare the Stieltjes transforms corresponding to X k and X k _\ with a third one, corre- 
sponding to the matrix Q hereafter with vanishing fc-th entry: noting k = aN + b (a G [0, TV — 1], b € [1, NJ) 
we define the following N x N matrices: 

v = v ab e ab = (a, b)e ab , (3.49) 

u = u ab e ab = X (2 \a, b)e ab , (3.50) 

Q = X k _ 1 -v = X k -u, (3.51) 

R={Q*Q-wiy 1 (3.52) 

ft = (QQ* - w/)- 1 

5 = {Xl^Xk.x - wiy 1 (3.53) 

T = (X* k X k - wl)- 1 (3.54) 



Then (3.481 holds if we can prove that 

Tr R - E Vab Tr S = F + (log N) c OiN" 2 ^ 1 ), (3.55) 
Tri? - E Uab TrT = F + {logNf OG/V" 2 ^" 1 ). 

holds with £-high probability where E t , ab (resp. E Uo6 ) means an integration only with respect to the random 
variable v ab (resp. u ab ), and F is a random variable, identical for both equations. We prove the first equation, 
the proof of the other one being obviously the same. For this, we want to compare S and R. Denoting 

U ={l + RQ*v)~ 1 R(l + v*QR)-\ 

we first compare Tr U and Tr R, and then Tr S and Tr U. For this first comparison, we introducing the 
notations 

P = {RQ*) ba Vab, q = (QR)abvl a , 

appearing in the following identity obtained by expansion^] 
Tr U = Tr(l + v*QR)- 1 (l + RQ*v)~ 1 R = ^ (-l) k+l Tr ((v*QR) k (RQ*v) l R) 

= Tr R + J](-l) fe Tr ((v*QR) k R) + ^(-1)' Tr ({RQ*v) l R) + £ (-l) fe +* Tr ((v*QR) k {RQ*v) l R) 

fc^l 1^1 k,l^l 

TrU = TrR V b^ QR ^ ab ( R2 Q*)baVab (QR 2 Q*)aa\Vab\ 2 Rbb , . 

1 + q l+p (l+p)(l + <?) 

For this last equality, we extensively used that if ^ then — (a, b). 
To compare now Tr S and Tr U, we introduce the notation 

fl = v*(QRQ* - l)v. 

A routine calculation yields X'^_ 1 X k -i — w = U^ 1 — f2, hence 

S = {l-Uny 1 !/. (3.57) 



3 All the expansions considered here converge with high probability. Anyways, they aim at proving identities between rational 
functions, which just need to be checked for small values of the perturbation. 
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As QRQ* — I = wlZ, we have fLj = rdibSjb where 



r = \v ab \ 2 wTl aa . 



Consequently, expanding (3.571 we get 

Sii = U-, 



so after summation over i = j, 



, -r U ib rU bj + U ib rU bb rU bj + U ib rU bb rU bb rU bj 
Tr S = Tr U + [(7 2 ] b6 r(l - CW)" 1 



(3.58) 



An argument similar to the one leading to (3.561 yields that, for any matrix M, 

M bb 



[(1 + RQ*v)- x M{\ + v*QR)~ v 



which yields, in our context, 
Rbb r 



U bb 



r(l +}>)(! + g) 



(l+p)(l + <?)' (l+p)(l + g)-ri? bh 



(l+p)(l + g)' 

(i?(l + v*QR)- 1 (l + RQ*v)- l R) bb 



(l+p)(l + 9) 



The numerator of this last expression is also, by the same rcasonning leading to (3.56) 

(R 2 Q*) ba V a bRbb (QR 2 Q* 



(R 2 )bb 



RbbVl a (QR )ab 



aa\Vab\ 2 Rl b 



1 



1+P 



(l+p)(l + g) 



Substituting the above expressions in (3.581 and combining it with (3.561, we get 

)bb 



Tr S - Tr R 



r{R 2 ) 
(1 + p)(l + q) - rR bb 

(l+p)(l + g) 
(l+p)(l + g) -rR bb 



vL(QR 2 )g b 

1 + 9 



(R 2 Q*)baVab , (<3-R 2 Q*)oa|« fe| 2 -Ri 



/;6 



1+P 



(l+p)(l + «) 



(3.59) 



In the above formula, only appears through p, q (in a linear way) and r (in a quadratic way) . All 
other terms can be bounded thanks to the estimates for S from Lemmas |3.6| and |3.7| which also holds for R, 
as stated in Remark |3.8| More precisely, the following bounds hold with C - high probability, for any choice 



of a and b (including possibly a = 6), and under the assumption (2.191: 



\(RQ*)ba\, \(QR)ab\ 




(log 


;N) C 


by ( 


3.38 


\Rbb\, \Tt aa \ 




(log 


;N) c \w\-^ 2 


by ( 


3.19 


\(R 2 )bb\ 

\{QR 2 )ab\ 
\(QR 2 Q*)aa\ 




(log 


; 7V)C| w |-i/2 ?? -i 


by ( 


3.40 




(log 




by ( 


3.41 




(log 


; 7V)C| w |i/2 ?? -i 


by ( 


3.42 



(3.60) 



Therefore, if we make an expansion of (3.59) with respect to v ab , we get for example for the first term of 
the sum in (3.591 with ^-high probability (remember that v ab satisfies the subexponential decay property) 



r(R 2 ) 



bb 



= / + O (u>CR 2 )66ft«„M 4 max {(RQ*)la, (QR)lb, ™RbbTl aa )) = f + O 



N 2 r] 



(l+p)(l + q)-rR bb 

where / is a polynomial of degree 3 in v a b- A calculation shows that the same estimate, of type g + O ^^2^ 
with g of degree 3 in v a b, holds when expanding the second term of the sum (3.591. This finishes the proof 



of (3.55). 



□ 
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Proof of Lemma \2.1S\ First since Theorem 1.2 holds in the Gaussian case (cf. Theorem 4.4 1 , then with 
(2.221, the estimate on the smallest eigenvalues (Lemma 2.9 1 and the largest eigenvalue (Lemma 2.10), the 



estimate (2.30) holds in the case of centered and reduced complex Gaussian entries. Furthermore, using 



(2.321, we deduce that Lemma 2.13 holds in the Gaussian case. We will therefore use the same method as in 



the proof of Lemma |3.9| replacing the matrix elements one by one to extrapolate from the Ginibre ensemble 
to the general setting. We assu me tha t = Xq is the Gaussian case and = X N 2 is the ensemble for 
which we want to prove Lemma 



2.13 



We know that Z^ 1} -< N Ve . and we even know that, for any fixed 



p > 0, Ed-Zj/^) c | p ) ^ N Cep (Theorem 4.4 is proved by bounding the moments). We will prove that, for any 
fixed p £ 2N and N sufficiently large, 



E (Z 



■(f) \ 

X K 2,C> 



< C p N ds v, 



(3.61) 



for some C > C independent of N and p. This is sufficient for our purpose, by the Markov inequality. 



From Xk and Xk-i, we defined the matrices R, S, T, in (3.521, (3.53), (3.541, From equation (3.59), one 
can write 

TrS* - Tri? = P w , z< Q(v ab ) + B», !2 ,s^ t , 

where P w . z .q is a degree 3 polyn omial whose coefficients depend only on Q. Moreover, a careful analysis of 
(3.591, using the estimates (3.60), proves that the coefficients of P are 0^(?/ _1 ), and that B is also Ox(f7 -1 )- 
0^(N~ 2 ), hence we obtain 



Note that v* b 



7 (f) _ 7 (f) _ 



a/(0 / xinW{E) (^N-^dEdnd^ 



+ J A/(0 J iX (v)<t>'(E)Re(P w , z , Q (v ab ) - P WtZtQ (0))dEdr]d(dZ 
Noting that |<//(-E)| ^ (1 + log E)E~ 1 \ i \ +>E>N z-2 > the first term in the above sum is 



N- 2 / |A/(0| / ^E^dEdrjd^ = {N- 2 \Af\ L1 ) = (N' 2 ) , 

y J J /n{4A + >£>A^- 2 } ) 

where we omit the dependence in / in the previous and next estimates, as / does not depend on N. 
Concerning the second term, it is of type Re"P f,Q(v ab ), where V has degree 3, vanishes at 0, with coefficients 
of order 0^(1) being independent of v ab and u ab . We therefore have 



7 (f) 

y(f) 



7(f) 
J Q,c 

7(f) 



A fc _i, A fe _! =V fiQ (v ab ) + 0^(N- 



(3.62) 



Z%> c + A fe , A k = P f , Q (u ab ) + 0^(N- 2 ) 



We can decompose 



(4L, c ) p - (4 / fc ) , c ) p = E Q (4'c) j '(Ari - Ar )■ 



3=0 

Since V has no constant term and the first three moments of v ab and u ab coincide, we get 

p-i 



E | 



((Z { £_ uc r) - E ((4^c) P ) | < E E |(4^c) J '| 0^(N- 2 ) < TV- 2 (O^(l) +E ((Zg)')) , (3.63) 



J'=Q 
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where we used that p is even, E 



(z ( f ly 



E 



< E 
j/p 



(z ( n f) c r 



j/p 



N <E 



and 



(3.64) 



Moreover, from (3.62), 



e (zgr a (4i 1>c ) p +E 



p 



E 



^((41 1>c ) p )+E 1 E ( 

E ((^ 3 c) P )<(04l)+E((4i 1 ,c) P )): 



7 (/) 



4 /} 



|A fc -i| J 



E(|A 



fc-i 



where we used Holder's inequality and the trivial bound A^-i = 0^(1) in the last equation. Using the last 
bound and (3.631, we obtain for any 5 > 

e ((4L,cr) - e ((42,c) p ) | < ^ 2 (n s + E ((41 I!C ) P )) ■ 



It implies 



E ((Z ( £j p ) + N s < (1 + iV- 2 ) (E ((4L,c) P ) + ^) 



As E ( (i?^/ c ) p ) ./V Cep and we obtain ( 3.61 1 by iterating ( 3.65 1 . It completes the proof of Lemma 



2.13 



(3.65) 
□ 



Proof of Lemma \2. 1J\ Notice in the previous proof, the third moment condition was first used in (3.63l. 
So all equations up to and including (3.621 are still valid. Introducing the notation Y = Z/y/N, we have, 



instead of (|3.63|, the following bound 
E 



p-i 



((41Lc) p ) -E((4{;y)| <5>|(y^|o^(Ar 2 ) ^iv- 2 (o x (i)+E((y£^)) . (3.66) 



3=0 



Following the rest of the argument in the proof of Lemma 2.13 we have proved that 



and this completes the proof of Lemma |2.14| 



(3.67) 

□ 



4 The local circular law for the Ginibre ensemble 

In this section, we derive the local circular law on the edge for the Ginibre ensemble ( |1.1[ ). This is required 
in the proof of Lemma |2.13| which proceeds by comparison with the Gaussian case. Along the proof, we will 
need the following partition of C, with distinct asymptotics of the correlation function Kn for each domain 
(in the following we will always take en — 

fil = {\A < 1 - e N }, the bulk, 

il 2 = {\z\ > l + e N }, 

fi 3 = C - (fii U Q 2 ), the edge. 
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We first consider the case when the test function is supported in f2i (this step is required even for the 
circular law on the edge, as the support of the test functions overlaps the bulk), then when it can overlap 
fi 3 . The first case is directly adapted from (lj, the second one requires some more work. 

In the following, we note cn(z) = X^fco fr ^ or ^ ne partial sums of the exponential function. 

4.1 The bulk case. We prove the following local cirular law for the Ginibre ensemble, with some more 
precision: the local convergence towards the Gaussian free field holds in the bulk, generalizing the global 
convergence result obtained in [23] , 

Theorem 4.1. Let < a < 1/2, f be a smooth non-negative function with compact support, and Zq be such 
that f Zo is supported in Q±. Let fix, . . . , fj,^ be distributed as the eigenvalues of a Ginibre random matrix 
(1.1). Then the local circular law (1.6) holds. 

More generally, if f additionally may depend on N , such that ||/||oo *S C, ||/'||oo ^ N , the local 
circular law holds when at distance at least ejy from the unit circle in the following sense. Define f^ N ^(z) — 
f(N a (z—Zo)), for Zq (depending on N ) such that f( N > is supported on flj. Denote cumjv(£) the l-th cumulant 

ofX [ } N) = E/ W (w) - N 1 ' 2a ^J c f(z)dzdz, and denotF\s N = TV 3 ^ Vol(supp(V/W))||/||°o ■ Then as 
N -+ oo, for t > 1, 

cumv(/) = / ^/|V/(z)| 2 dA(z) + 0(M if£ = 2, 

carnal) | Q ^ if £^2. [ > 

In particular, noting a 2 — J |V/(z)| 2 dA(;z), if Sn *C a 2 , the linear statistics ijfy converge in law to a 
reduced centered Gaussian random variable. 

As a first step in the proof, we need the following elementary estimate. 

Lemma 4.2. Let 

k N ( Zl ,z 2 ) = ^e-T (W+W-a.iSi). (4 .2) 
There is some c > such that uniformly in \zi~z^\ < 1 — we have 

K N ( Zl ,z 2 ) = k N (z u *2) + 0(e- c ^ N ^), 



where Kjq is the Ginibre kernel (1.2). 

Proof. This is elementary from the following calculation, for any \z\z 2 \ < 1: 



a 

£>N 



= k N (z u z 2 ) + [iVe-f(l-l 2 +l-l 3 ^^ £ \*M 
= k N {z u z 2 ) + o(Ne' N ^ + \^' 2 - 1 -^ z ^ l — , 

V \-\Z\Z-z\ i 

K N ( Zl ,z 2 ) = k N ( Zl ,z 2 ) + fiVe-T(i-l^l) 2 1 ) (4.3) 

V l-\z lZ2 \J 



MM™ :=max(||/'||| 0) ||/'|| o||/"||ao,||/W||oc). 
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□ 



where we used Stirling's formula and the fact that (\zi\ 2 + | ^2 1 2 ) / 2 — 1 — log \ziZ2\ > \z\z 2 \ — 1 — log \z\Z2\ > 
\{l-\ziz 2 \f. 

The next following estimate about Kjq will be used to bound it for any distant Z\ and z%. 
Lemma 4.3. There is a constant c > such that for any N G N* and any Z\, z 2 in C, 



,-JVi 



\Z\Z% 



AT+1 



1 



N\l 



\K N { Zll z 2 )\ < cN 

Proof. In the case \z\Z2 — 1| < 1/V-/V, by case (i) in Lemma 

N 

K N (zi 7 z 2 ) = —e 



-N 



i*ii z +i*2r 



4.9 



Zl^2| 

we get 



e Nz lZ2 I _ erfc (ViV^(z 1 ^)) + 0(7V- 1/2 ) 



As erfc is uniformly bounded in C, we get \K^{zi, z 2 )\ = 0(e 



— N- 



)• 



In the case -7= ^ l^i^a — 1| < S, for 5 fixed and small enough, the estimate was obtained in pi (cf. the 
proof of Lemma 8.10 there, based on Lemma 4.9 here). Finally, in the case \z1_Z2 — 1| ^ 6, an elementary 
calculation from Lemma [4. 1 1 1 implies that the result holds, no matter that Z1Z2 is in D — U or D c — U 



□ 



to imply the weak convergence to the Gaussian distribution, and it also yields the local circular law (1.6): if 



Proof of Theorem 4-1 Proving equation (4.1) is sufficient: this convergence of the cumulants is well-known 
to imply the weak convergence to the < 
tojv(£) denotes the £-th moment of Xj 

m N (£) = 



one can write 



E 



A, 



cumjv(ii) . . . cumjv(^) 



for some universal constants A, 1) ... 1 ^'s, where the ij's are in N*. Hence writing $ = [a 2 



m N (£) = O 



Sn) 1 ^ 2 we get 



Consequently, for any e > 0, 



P(\Xf> \ > N e <f>) < m N {l)N- 



0(N~ t£ ) 



which conclu des t he proof of the local circu lar law by choosing £ large enough 



To prove (4.1l, first note that from 



4.2 



evaluated on the diagonal z\ = z 2 , cuhiat(I) = 0(e~ c ( logAr ' ), so 
we consider now cumulants of order £ J? 2. Due to a peculiar integrability structure of determinantal point 
processes, the cumulants have an easy form for the Gaussian matrix ensembles, as observed first by Costin and 
Lebowitz 10 , and used later by Soshnikov |27|. For our purpose, in the context of the Ginibre ensemble, we 
will use the following very useful expression (4.5 1 due to Ameur, Hedenmalm and Makarov |;1| , which requires 
first the following notations. The usual differential operators are noted d = \{d x — id y ), 8 = \(d x + id v ), Ag = 
didi + ■ ■ ■ + dgdi, we will also make use of the length-^ vectors zlg = (z, . . . , z), h — (hi, . . . , he), and the 
following function often appearing in the combinatorics of cumulants for determinantal point processes: 

£\ j 



.z e ) 



3=1 J fciH hfei=«,fei, 



ki\. 



.kA — 
J m=l 



n f [N) ^ h 



Then the £-th cumulant of X^ is then 



cumjv(^) = / F t (zi, zi)K N (z 1 ,z 2 )K N (z 2 , z 3 ) . . . K N (z t , z x ) dA(zi, . . . , z e ), 



(4.4) 
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Remarkably, some combinatorics can prove that Fg(ztg) = 0: it vanishes on the diagonal. Moreover, the 
following approximations of Ft close to the diagonal will be useful: 

F i (zl e + h) = Y / T j (z,h), where T j (z,h) = V (d a B^F e )(zt e ) — — . 

We will also need the third error term after second order approximation, r(A, h) = Fg(z1e + h) — T\(z, h) — 
T 2 (z, h). The last notation we will need from [l is Zg(z) — Y^i<j (pSjFi) { z ^-l)i an d the area measure in C e 
will be noted dA(z, hi, ... , hi) = dA(z) dA(ft-i) . . . dA(ht), where dA(z) = d 2 z. Then, as proved in [l], for 

£ ^ 2 

cum N (£) = A N {£) + B N (£) + C N {£) + D N {£) + E N {£) (4.5) 

where 

A N {£) = [ ft [YXd i d j F £ )(zl £ )h 1 h 2 \ K N (z,z + hi)K N (z + h x ,z + h 2 )K N (z + h 2 ,z) dA(z,h 1) h 2 ) 

B N {£) = ^ [ y2(d t 2 F e )(zt e )h 2 1 K N (z,z + h 1 )K N (z + h 1 ,z)dA(z,h 1 ) 

C N {£) = 2 (Zi(z)hih 2 ) K N (z, z + h 1 )K N (z + h x ,z + h 2 )K N (z + h 2 , z) dA{z, hi, h 2 ) 



D N (£)= / (A e F e )(zl e )\hi\ 2 K N {z,z + hi)K N {z + h 1 ,z)dA(z,h 1 
Jc 2 



En(£) = / r(z, h)K N (z, z + hi)K N (z + hi, z + h 2 ) . . . K N (z + hi, z) dA{z, hi,..., he). 

Remember that the support of / W is at distance En from the unit circle, so each one of the above integrands 
vanishes if z is out of the disk with radius 1 — £jy. Therefore, by using Lemma |4.3| to restrict the domain, 
and then Lemma |4.2| to approximate the kernel strictly inside the unit disk, one easily gets that 



cum N (£) = A N {£) + B N {£) + C N {£) + D N (£) + E N (£) + 0(e- c ^ N ^ 2 ), 
for some c > 0, where 

A N (£)= [ ^[Y J {d l d 3 F l ){zt,)h 1 h 2 \ k N (z, z + h x )k N (z + hi, z + h 2 )k N (z + h 2 , z) dA(z, hi, h 2 ) 



C3n{||/ t || 00 <e N } 



Bn(£)= & [ y2(d 2 F e )(zt e )hjk N {z,z + h 1 )k N {z + h 1 ,z)dA(z,h 1 ) 

JC 2 n{||/ l |U<e JV } i 

Cn{£) — 2 / 5i [Zg (z)hih 2 ) k N (z, z + hi)k N (z + hi, z + h 2 )k N (z + h 2 ,z) dA(z, hi,h 2 ) 

" / C 3 n{||/ l || 00 <e JV } 

Dn{£)= I (A e F e )(zt e )\h 1 \ 2 k N (z,z + h 1 )k N {z + h 1 ,z)dA(z,h 1 ) 

JC 2 n{||/7.|U<e JV } 

En{£) = I r(z, h)k N (z, z + hi)k N (z + hi, z + h 2 ) . . . k N (z + h e , z) dA(z, hi,..., hi). 

i y c c+in{||) l || 00 < ejv } 



2G 



Following closely the proof in [l], integrating the hi's for fixed z, the terms An(£), Bn(£),Cn(£) are 
O ( e -c(io g A0 2 ^ Moreover, [l] proved that, for £ ^ 3, AgF( vanishes on the diagonal, and that A 2 F2(z l7 z 2 ) | Zl=z . 
i|V/W(»| 2 . Consequently, Dn(£) = if O 3 and 

^(2) = s / \Vf {N Hz)\ 2 \h\ 2 ^e- N W 2 dA(z,h) + 0(e-^^ 2 ) 



2 



i / |V/W(z)| 2 dA(z) + 0(e- c ( lo ^) 2 ) = -^||V/||2 + 0(e- c ( 10 ^) 2 ). 

47T ./r 47T 



To finish the proof, we need to bound the error term E^(£). We divide the set of possible (z,h) in two 
parts. Following 111, define Y n g — {(z,h) £ C e+1 : z £ supp(V /^)}. On the complement of Y n ^, 
z £■ supp(V/W), and then r(z,h) = 0. On Y n/ , r(z,h) is of order 0(|| J^^ 3) |U H/iH^) = 0(^ 3 1|/l|Ue^), 
each kpf is O(N), and the domain of integration has size 0(Vol(su PP (V/W)) 

E N {£) = O ((log NfN 3a ~i\l fl^ Vol(supp(V/W))) , 
concluding the proof. □ 

4.2 The edge case. In the edge case, i.e. |zo| £ O3, we have the following theorem for the Ginibre ensemble. 

Theorem 4.4. Let /ii, . . . , ^jv be distributed as the eigenvalues of a Ginibre random matrix and zq £ Q3. 

Suppose that f is smooth and compactly supported, and let f Zo {z) = N 2a f(N a (z — z )). Then for any 
< a < 1/2, the estimate [1.1) holds. 

We begin with the proper bound on the first cumulant, noting as previously f^ N \z) = f(N a (z — zq)). 
and Xf> = E/(^ - zo)) - f f D f^(z) dA(z). 

Lemma 4.5. With the previous notations, for any e > 0, E(x| ) = 0(N~^ +£ ) as N —} 00. 
Proof. From the definition of the 1-point correlation function, 

From Lemma |4.2[ as / is bounded, 

N[ f( N \z)(j j K N (z,z)-l)dzdz 



dzdi 



!2i 



= N /W(z) (jjk N (z,z) - dzdz + O (e-< l °^ = O (e^ 1 ^) 2 ) 



where kjy is defined in (4.2). On f2 2 , by Lemma 4.3 K^{z,z) — 0(e c ^°s N ) ), Consequently, 



fW(z)K N (z,z)dzdz = O (e- c < lo s") ! ) 



Without loss of generality, we consider the case \zq\ = 1 for simplicity of notations. On we now use 
Lemma |4. 10| which yields that, uniformly on — Sn < £ < £n, 

N 

7 V~ £ " u ' V^i 



K,,fi - ! !- -:) = -( F n f ^ erfc(ViV - l M (f)) + ) ) . (4.(5) 
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where t = + e) 2 - Take any smooth function of type g(l +e) = G(N a e). By polar integration, we just 

need to prove that the following quantity is 0(N a ~i +e ): 

J " (1 + e)g(l + e) \K N (1 + e, 1 + e) - de 

(1 + e)(g(l + e) - «?(!)) (k n (1 + e, 1 + e) - ^\ de + O ( e -<=0°s"> 2 ) , (4.7) 



where in the last equality we used the simple yet useful facts that J c Kn(z, z) dA(z) = 1 and \Kn(z, z) — 

f l| z |<i| = O(e" c ( lo s jv ) 2 ) uniformly out of fl 3 . We now can write g(l + e) - g(l) = N a G'(0)e + 0(N 2a e 2 ), 

and use (4.6) to approximate Kn(1 + e, 1 + e) — Using the fact that fi(z) = + 0((z — l) 2 ) for z 

close to 1, still denoting t = + e) 2 , we have fi(t) — \/2e + 0(e 2 + A -1 ), so the integral term in (4.7) 

can be written 



A 
n 



(1 + O(e))(A a G'(0)e + 0(A 2a e 2 )) Q + 0(e + A" 1 )^ 

erfc(%/2A|e| + 0(A~ 1/2 + A 1/2 e 2 )) (l + 0(A~ 1/2 )) de 

One can easily check by bounding with absolute values that all terms involving O's , except the one in the 
erfc function, contribute to 0(A a ~5 +e ), so we just need to prove that 



A / eerfc(v / 2]V|£|+0(iV- 1/2 + A 1/2 e 2 ))de = 0(iV-5). 



As erfc' = 0(1) uniformly, this is the same as proving A Jf ^ e erfc(v / 2A|e|)de = 0(A 1 ^ 2+E ), which is 
obvious as it vanishes by parity. □ 

We now prove the convergence of the second moment. 

Lemma 4.6. Under the hypothesis of Theorem\4-4\ as N ^ oo, 



cum w (2) = i- J |V/(z)| 2 l^dA(z) + ill/Til^ +o(l) 

where we use the notation A = {z | Ji(zzjj) < 0}, /t( z ) = f{i- z o z ) and WgW^/2 = 4^2 J R 2 f ^^-^^ ) daidy. 
Proof. Note that by polarization of formula (4.4), 



cuTn N {2) = \ f (/W(z 1 )-/W(z 2 )) 2 |A w (z 1 ,z 2 )| 2 dA(z 1 ,z 2 ). 



One easily gets 



cum w (2) = i / (/<">(*i) - f^(z 2 )) 2 \K N (z u z 2 )\ 2 dA< :, . 



+ \ f (f W M f iN) {z 2 )) 2 \K N { Zll z 2 )\ 2 dA(z 1; z 2 ) + O (e-^) 2 ) 
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the contributions of the domains C x fl 2 and f^i x being easily bounded by Lemma 4.3 For the Q 2 term, 
one can easily reproduce the method employed in the bulk, in the previous subsection, approximating Kn 
by fcjv thanks to Lemma 14.21 to obtain 



\ I (f^Hz 1 )-f^(z 2 )) 2 \K N (z 1 ,z 2 )\ 2 dA(z 1 ,z 2 ) = ^\\xA^f\\l + o(l) 

2 JQ.\ 47T 

The most tricky part consists in evaluating the £7§ term. To calculate it, we will need the following notations, 
where 1 = arg(zi), 9 2 = arg(z 2 ): 

(i) = \!^ {N Kz 1 )-f N \z 2 )) 2 \K N {z l ,z 2 )\ 2 dk{z 1 ,z 2 ), 

(ii) = 1 J n (f (N) (zi) ~ f (N) (z2)) 2 \K N (z u z 2 )\ 2 dA( Zl ,z 2 ), where SI = ({(z u z 2 ) € Of : \6 X - 6 2 \ > N'^}), 
where e > is fixed and will be chosen small enough, 

, z 2 



(hi) = i^(/W(e^)-/ (A,) (e i92 )) 2 |^(^i,^)| 2 dA(z 1 , 
(iv) = ^L ^'pliqr^^ -^'^^'-'^^kr^l-dA^,^), 

(vi) - JL f j/ilnl-JitfipV j 

^ _ 8tt 2 Jr 2 (r!-r 2 ) 2 ^1^2- 

We will prove successively that the six above terms differ by o(l) a N — > oo, which will conclude the proof 
as the last one is exactly IH/tIIhi/2- 

First, the difference between (i) because (ii) is o(l), because it can be bounded by 



ejvJVII/USo sup / 1 \z 1 -z 2 \'\K N (z 1 ,z 2 )\ 2 dA(z 2 ), 

Zi6fi 3 J\z 2 -z 1 \<N~2 +e 

where N~ a eN corresponds to the position of z\ such that \6\ — 2 \ < N~^ +£ and \f^(zi) — f^ N H z 2 )\ ^ 0, 
and iV 2a |]/ , ||^ comes out by using the Lipschitz property for f( N \ Moreover, From the lemmas 4.9 and 
4.11 it is easy to prove that uniformly in C 2 



so finally one can bound the difference between (i) and (ii) by 

^AHI/IIL / . r 3 ( N ) dr = O {s^N^WfWl) = o(l) 
Jo<r<N-2+ E \l + VNrJ 

when e < | ( ~ — aj . 
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Now, the difference between (ii) and (iii) is o(l) as well: by an order 1 approximation around e 1 ® 1 and 
e 1 ® 2 , this difference is of order at most 

£Jv||/ (JV) '|U|/ (iV) ||oo / \K N (z 1 ,z 2 )\ 2 dA(z 1 ,z 2 ) 
^n({/(»)( Zl )/o}u{/C)(z 2 )|/o}) 

r n 2 

^e N N a / — dA(z 1; z 2 ) 

^nn({/W(*i)^ }u{/W(%)|/o}) (l + VN\e 10 ^ -e lS =|) 2 



dr 



r>N -i/2+e (1 + y/Nr) 



^e 6 N N-2~ s , 

where the derivatives are radial ones, so (ii) — (iii)= o(l). 

For the difference between (iii) and (iv), we want to get a good approximation for Kn{z\, z<i). Note that, 
as 1 6*i — 6*2 1 3> N~ x / 2+e 3> 2en, one easily has that on Q, | arg(ziZ2 — 1)| ~ ? uniformly, in particular the 



formula (ii) in Lemma 4.9 provides a good approximation for e^-iiN z\z 2 ): together with the asymptotics 



erfc(z) ~^ f^r, for | argz| < *f , and noting that fi'(z) = 2 ^ ( g Zlf2 , we get 



ejv_ 1(^1^2) ~ , — — — -e N (z 1 z 2 ) N , 

N-yoo v^riV( Zl z2 - 1) 



so in particular 



l^v(^ 2 )| 2 ~ - 2 , 1 - 1 1|a 9 e- N ^ 2 ^% lZ2 \ 2N 

JV->oo 7T Z 27T|ZlZ2 — 1| — 2 

_ ^ 1 -JV(|* 1 | a +|* a | a -2)i .2tf 

w^co 7T 2 27r|e ie i - e ie H 2 ' 

the last equivalence relying still on the fact that, on O, \0\ — 2 \ 3> |zi — e' 6 * 1 1, |2T 2 — e l6 * 2 1 . This proves that 
(iv)-(iii)=o((iii)). 

To obtain (v)-(iv)=o((iv)), just note that \ Zl \ 2 - 1 - log(|zx| 2 ) = \{\ Zl \ 2 - l) 2 + 0(||^i | 2 - 1| 3 ), so by a 
simple saddle point expansion we get 

xe- Nx2+N x 2N dx ~ / xe-^^-D'dx ~ - / e^du = " 



By integrating (iv) along the radial parts of zi,z 2 , we therefore get (v), up to the domain where \9% — 2 \ < 
N- which can be shown to have a negligible contribution, as previously. 

Finally, (v) converges to (vi), as a simple consequence of a change of variables and a dominated conver- 
gence. □ 

For all cumulants of order I ^ 3, we begin with the following very crude bound. 

Lemma 4.7. Under the conditions of Theorem D.4-, for any fixed function f, e > 0, and I 3, we have (as 
N -> 00) 

cum N (£) = 0(N 3 ). 
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Proof. In the decomposition (4.5 1, it is obvious that the terms An(£), Bn(£),Cn(£), Dn(£) are 0(7V 3 ), just 



by bounding Kn by N in fii U and using lemma |4T3| if one point of Kn is in f^- Concerning the term 



En(£), by reproducing the argument in the proof of Theorem 4.1 we just need to prove that it is 0(N 3 ) 



when all points in the integrand are restricted to O3, i.e. it would be sufficient to prove that 

\K N (z , Z!)K N {z, zi) . . . K N (zi, z )\ dA(z , ...,z t ) = 0(A 3 ). 



-+1 



using the inequality (4.8 1 and integrating along the width of W3, we just need to prove that for some e > 0, 



N 1 -^ [ — ... -J dx Q ...dx t = 0(N 3 - £ ). (4.9) 

J l + ^N\ XQ - Xl \ 1 + ^N\ Xe - XQ \ 

A simple calculation yields, for any — 1 < o < 6 < 1, 

-du — O 



(-1,1) 1 + VN\a -u\l + VN\u -b\ V VN J 1 + y/N\a - b\ ' 



so by integrating successively the variables X2, ■ ■ ■ ,xe in (4.9 1, we therefore obtain that a sufficient condition 

18 N /(-l,l) a (l+VN\ 1 xo-x 1 \) 2 dX ° dXl = °( 7V3_£ )' Which 18 obviouS - D 

Relying on the Marcinkiewicz theorem, the above initial bounds on the cumulants can be widely improved. 

Lemma 4.8. Under the conditions of Theorem D.4, for any fixed function f,e>0, and £ 3, we have (as 
N -> 00) 

cum N (£) = 0{N e ). 

Proof. Let rW = N~ e xi N \ Let cum N (£) be the £-th cumulant of 7^: cumjv(/) = AT" fe ciim^f). 



Hence, as a consequence of Lemma 4.7 there is a rank such that for any £ > £q, cunijv(^) — > as N — >• 00. 
We wish to prove that 

sup cum w (f) — > (4.10) 

as well. Assume the contrary. We therefore can find 6 > independent of N, an index j\ 6 [3, £0] and a 
subsequence (AT*.) such that 

\cwi Nk {jx)\ > 8, Icumjvjji)! 1 ^ 1 = sup {IcumjvJ^I 1 ^}. 

Consider now the random variable = Y^ Nk ' /|cvlmjv fe 0'i)| 1/jl - Then all the cumulants of go to 
except the Ji-th one, equal to 1, and eventually some other cumulants with indexes in [3,4)], which are all 
uniformly bounded. Let /j,k be the distribution of Z^ k \ As the first and second cumulants of = are 
uniformly bounded, K(Z^ ) is uniformly bounded so the sequence (/Xfc)fcj>o is tight. 

We therefore can choose a subsequence of (/ifc), called (ni k ), which converges weakly to a measure v, 
and we can moreover assume that all of its cumulants (which are uniformly bounded) with indexes in [3, £qJ 
converge. Let Z be a random variable with distribution v. As fj,i k has any given cumulant uniformly bounded, 
it has any given moment uniformly bounded, so by the Corollary of Theorem 25.12 in 15], Z has moments of 
all orders, and those of Z 1 - 1 ^ converge to those of Z. Hence Z has cumulants of all orders and those of Z^ %k ' 
converge to those of Z. Hence v has all of its cumulants equal to for £ > £ , and its cumulant of order 
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ji 6 [3,^o] equal to 1. From the Marcinkiewicz theorem (cf. the Corollary to Theorem 7.3.2 in [20]), such a 
distribution v does not exist, a contradiction. 

Equation (4.101 therefore holds, so we proved that for any arbitrary e > and any £ 1, cuni/v(^) = 
0(N £e ). For fixed £, by choosing e small enough, we get the result. □ 

Finally, we note that the previous bounds on all cumulants of X^, each of order at most N £ , allow 
to conclude the proof of Theorem |4.4| by the Markov inequality. The explicit asymptotic form of the 
second cumulant (Lemma 4.6 1 is more than what we need, we made the calculus explicit as one expects that 
Xj N ^ converges to a Gaussian centered random variable with asymptotic variance J |V/(z)| 2 l z£j 4dA(2:) + 
sII/tIIhi/s- 



4.3 Estimates on the partial exponential function. For the following lemma we need the function erfc(z) 
2 



3= f +00 e^dw. 



Lemma 4.9 (Bleher, Mallison |6|). For S > small enough, let (J,(z) = yz — log z — 1, be uniquely defined as 
analytic in D(l, S) and + x) > for < x < S. Then for any M > 1, as N — > 00, e~ Nz eN{Nz) has the 
following asymptotics: 



2?r . 

3 > 



(1) \ erfc(^ M (z)) + 0(N-^) if \z-\\< ^L; 

M ubiz) erfc (^M(^)) (l + O ((^pv)) if 7= < \z - IK § and | arg(z - 1)| < 

^ 1 " 2vfe) erfc(-^M(^)) (l + O (j^)) if ^ ^ \z - 1\ ^ 6 and \ arg(z - 1) - tt| < f . 

In case of real z, the above estimates can be gathered as follows. 
Lemma 4.10 (Wimp |9|). Uniformly for t ^ 0, 



ejv(JVt) = l <t<i + ^ ^ erfc(ViV/x(t)) ^1 + O _ 



where fi(t) = y/t — log t — 1 is defined to be positive for all t, contrary to Lemma 4. 9 

Note that in the above lemma, the asymptotics are coherent at t = 1 because 

1 fi(t)t f -| as i -> 1" 
1 



y^2 t - 1 [ ^ as(4l+ ' 

Lemma 4.11 (Kriecherbauer, Kuijlaars, McLaughlin, Miller [19"). For < a < 1/2, let U = {\z — l\ < N~ a }. 
Then there exists polynomials hj of degree 2j such that for all r € N* we have 

f e Nz (1 - e Ar z Ar e- W2 F JV (z)) for z e D - U, 

[ e z Fjv(z) for z E D — U, 

where 

1 / ^ fy(z) 



" ^(1 - z) 1 1 + g W - + ° V A"- I J - :!■- 



uniformly in C — U . 
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